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摘要 : 本 文采 用 热 局 部 非 平衡 模型 对 完全 填充 金属 泡沫 的 水 平 圆 管 内 的 流 场 、 温 度 场 进行 了 数值 模拟 , 研究 了 金属 泡沫 的 关键 
结构 参数 如 孔径 d, 和 孔隙 率 s 等 对 管内 流动 与 换 热 特性 的 影响 。 研 究 结 果 表明 ， 填 充 金 属 泡沫 后 ， 管 内 流 场 变 得 比较 均匀 ， 金 
属 泡沫 圆 管 壁面 处 的 速度 梯度 远 高 于 光 管 ， 壁 面 边界 层 明 显 变 薄 ;孔径 对 压 降 的 影响 大 于 孔隙 率 的 影响 ， 当 孔径 越 大 ， 和 孔隙 率 
越 高 时 ， 压 降 越 小 ; 反之 ， 当 孔径 越 小 ， 孔 隙 率 越 低 ， 换 热 越 强 ， 综 合 考虑 金属 泡沫 参数 对 换 热 和 压 降 的 影响 ， 在 本 文 的 研究 
工 况 下 ， 最 佳 的 孔径 范围 是 lmm-2mm。 
关键 词 : 金属 泡沫 ， 完 全 填充 ， 热 局 部 非 平衡 
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A Local Thermal Non-Equilibrium Analysis of Hydrodynamic and Heat 


Transfer Characteristics in a Tube Filledwith Metal Foam 


WANG HuiGUO Lie-Jin 
(State Key Laboratory of Multiphase Flow in Power Engineering, Xian Jiaotong University, Xian, China, 710049) 

Abstract: A local thermal non-equilibrium model has been considered for the case of a tube fully filled with metal 
foam. The influence of metal foam key parameterson hydrodynamic and heat transfer characteristics has been 
investigated, such as pore size and porosity. The results show that the velocityin the metal foam tube becomes more 
uniform than that in an empty one, also a higher velocity gradient and thinner boundary layer are found. The pore 
size has a moreobvious influence on the pressure drop than porosity does. When the metal foam has a bigger pore 
size and a higher porosity, the pressure drop in the tube becomes lower. Otherwise, a smaller pore size and a lower 
porosity lead to a higher pressure drop, but an enhancement of the heat transfer. Considering the pressure drop and 
heat transfer both, the results indicate that the most appropriate pore size is 1mm-2mm. 
Key words: metal foam; fully filled; local thermal non-equilibrium 
0 前 音 并 考 虑 了 热 弥 散 效应 ， 获 得 了 与 实验 结果 比较 吻合 

的 模拟 结果 。Boomsma 上 等 以 水 为 冷却 介质 ， 对 压缩 
开 孔 金属 泡沫 由 于 具有 孔隙 率 高 、 密 度 小 、 比 。” 铝 泡沫 用 于 紧 竣 式 换 热 器 的 换 热 性 能 进行 了 实验 研 
表面 积 大 等 优点 ， 使 其 成 为 一 种 理想 的 强化 换 热 材 究 ， 指 出 在 相同 和 泵 功 下 压缩 铝 泡沫 换 热 器 的 热 阻 比 
料 ， 成 为 近 几 年 强化 换 热 的 研究 热点 。 最 好 的 商业 换 热 器 的 热 阻 低 2-3 倍 。Lu 和 Zhao5 
国内 外 已 经 有 很 多 学 者 对 金属 泡沫 强制 对 流 问 ”等 对 金属 泡沫 应 用 于 套 管 式 换 热 器 的 情况 进行 了 理 
题 进行 了 实验 研究 和 数值 模拟 。Zhao 和 Kimt"l 等 人 ” 论 分 析 和 解析 建 模 ， 针 对 内 管 与 环形 空间 进行 了 解 
在 真空 和 常 压 环境 下 ， 对 FeCrAIY 金属 泡沫 填充 的 析 研 究 ， 分 别 得 到 了 两 种 情况 下 的 解析 解 ， 通 过 解 
平板 通道 内 空气 的 强制 对 流 换 热 进行 了 实验 和 数值 。 析 解 得 到 的 数据 发 现金 属 泡沫 换 热管 能 大 大 强化 换 
模拟 分 析 ， 研 究 了 孔 聊 率 、 孔 密度 和 固体 导热 系数 。” 热 。 
对 换 热 的 影响 。Calmidi 和 Mahajanr ”对 铝 泡沫 填充 目前 针对 金属 泡沫 强化 换 热 的 实验 研究 多 集 
的 平板 通道 内 的 强制 对 流 换 热 进行 了 数值 模拟 和 实 。 于 平板 通道 内 的 情况 ， 而 数值 模拟 的 研究 中 能 量 方 
验 研究 ， 其 中 数值 模拟 过 程 采 用 热 局 部 非 平衡 模型 ” 程 多 采用 局 部 热平衡 模型 ， 即 认为 流体 相 的 温度 和 
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固体 相 的 温度 相等 。 这 种 假设 在 固体 相 的 导热 系数 
较 大 的 时 候 会 产生 比较 大 的 误差 。 本 文采 用 局 部 非 
热平衡 模型 对 完全 填充 金属 泡沫 水 平 管内 的 对 流 换 
热 进 行 数值 模拟 ， 研究 金属 泡沫 关键 结构 和 性 能 参 
数 对 换 热管 内 流 场 和 温度 场 的 影响 规律 。 


1 物理 模型 


本 文 所 描述 的 物理 模型 如 图 1 所 示 , 半径 为 2mm 
的 水 平 圆 管内 完全 填充 均匀 且 
平均 速度 为 wu 温度 为 7 的 热流 体 流 过 该 圆 管 及 
内 部 的 金属 泡沫 ， 冷 却 后 从 圆 管 的 右 端 流 出 。 圆 
壁面 具有 恒定 的 温度 7 且 Tvw<7Thn。 
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图 1 物理 模型 
Fig.1 Physical model 
2 数学 模型 


本 文 研究 中 采用 的 假设 如 下 : 

(1) 流体 不 可 压缩 ; 

(2) 忽略 自然 对 流 和 辐射 换 热 的 影响 ; 

(3) 忽略 热 弥散 和 粘性 扩散 的 影响 ; 

(4) 各 物性 参数 为 常数 。 

基于 以 上 假设 ,应 用 Darcy-Brinkman-Forcheimer 
动量 方程 和 局 部 非 热平衡 〈 | Thermal None 


Equilibrium ) 模型 ， 可 以 得 出 描述 该 模型 流动 和 传 热 
过 程 的 无 量 纲 控制 方程 : 
连续 性 方程 : 
0 1 0 
wj (]) 
z 方 向 动量 方程 : 
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向 同性 的 金属 泡沫 。 
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固体 相 能 量 方 程 : 
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为 使 得 计算 结果 具有 更 强 的 适应 
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FE， 定义 无 量 纲 参 数 如 下 : 
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其 中 ，Dea 为 达 西 数 ， 用 来 表征 多 孔 介 质 区 域 的 
渗透 能 力 。 渗 透 率 K 与 多 孔 介 质 的 结构 有 关 ， 可 根 
的 经 验 公式 求 出 ， 并 没有 固 
的 表达 式 。 本 研究 中 采用 Calmidi 吕 等 人 的 实验 拟 合 
公式 对 天 和 三 进行 求解 。 其 表达 式 如 下 : 
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式 中 , q, 为 金属 泡沫 的 平均 直径 ,df 为 泡沫 金属 内 金 


属 骨 架 的 纤维 丝 径 ， 它 与 四 的 关系 为 : 
d (1—é) 
En 1.18 (| (8) 


3 网 格 划分 


本 研究 所 用 网 格 如 图 2 所 示 。 在 流体 入 口 处 和 
近 壁 面 处 由 于 流体 的 温度 和 速度 梯度 变化 相对 较 大 ， 
计算 中 比较 容易 发 散 ， 为 较 好 的 捕捉 流体 的 参数 的 
变化 情况 ， 在 以 上 的 区 域 对 网 格 进行 加 密 处 理 。 同 
时 为 节省 计算 时 间 ， 其 他 区 域 采 用 相对 较 疏 的 网 格 。 


时 ， 从 壁面 到 管 中 心 的 固体 和 流体 的 温度 都 显著 增 
加 ， 但 固体 和 流体 间 的 温差 变化 不 明显 。 
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图 2 计算 区 域 网 格 划分 
Fig.2 Meshing method of computational domain 

数值 计算 中 ， 采 用 有 限 容积 法 对 所 建立 的 模型 
方程 进行 离散 求解 , 压力 和 速度 的 耦合 采用 SIMPLE 
算法 ， 对 流 项 采用 二 阶 迎风 格式 。 计 算 中 速度 和 温 
度 松弛 因子 设 为 0.7。 计 算 收 敛 的 判 据 为 系统 的 剩余 

质量 和 能 量 方程 的 残 差 分 别 降 到 le-8 和 le-10 以 下 。 

一 “4 计算 结果 及 讨论 

4.1 速度 分 布 
如 图 3 所 示 ， 为 不 同 孔 径 下 完全 填充 金属 泡沫 的 
水 平 管内 轴 向 速度 的 分 布 情况 ， 及 与 光 管 速度 分 布 的 
@ 对比。 由 图 3 可 知 ， 填 充 金属 泡沫 后 管内 速度 变 得 比 
较 均 匀 ， 金 属 泡沫 圆 管 壁 面 处 的 速度 梯度 变 大 ， 壁 面 


边界 层 明 显 变 薄 。 管 内 轴 向 速度 随 着 金属 泡沫 孔径 的 
增 大 而 增 大 。 这 是 由 于 在 孔隙 率 一 定时 当 金 属 泡沫 的 。 ”图 4 泡沫 参数 对 温度 场 的 影响 (a. 孔 径 ，b. 孔 阶 率 ) 
~ 孔径 增 大 时 ， 单 位 长 度 上 和 孔 的 个 数 减 小 ， 因 此 流体 受 Fig.4 Effects of foam structure parameters on 
a 到 的 阻力 相应 减 小 ， 此 时 管内 轴 向 速度 将 会 增 大 。 temperature distribution (a. pore size, b. porosity) 
10 pre 4.3 压 降 
>< 0 如 图 5 所 示 ， 为 金属 泡沫 结构 参数 对 管内 压 降 
Et i 4 影响 。 由 图 5(9) 可 知 , 当 孔 径 由 0.5mm 增 大 到 Imm 
一 | 时 ， 压 降 急 剧 减 小 ， 而 当 和 孔径 大 于 lmm 以 后 ， 压 降 
A 上 和 随 孔径 的 增加 继续 减 小 但 是 减 小 的 不 明显 。 由 图 5(b) 
让 交管 尾 ~ 可 知 ， 当 金属 泡沫 的 孔径 较 小 时 ， 随 着 孔隙 率 的 增 
| i 大 ， 压 降 和 急剧 降低 。 这 是 由 于 孔隙 率 增 大 时 阻碍 流 
i \ 体 流 动 的 固体 的 体积 分 数 呈 比例 减 小 ， 从 而 使 得 流 
M2 体 流动 的 阻力 减 小 。 但 当 孔 径 大 于 lmm 时 ， 和 孔隙 率 
UUn 的 变化 对 压 降 的 影响 不 明显 。 
图 3 不 同 孔径 下 金属 泡沫 填充 管内 的 速度 分 布 〈s=0.9) 35 
Fig.3 velocity distribution in metal foam tube under 0 a 
different pore size (s=0.9) 25| 


4.2 温度 分 布 | 
如 图 4 所 示 ， 为 金属 泡沫 结构 参数 对 填充 管内 的 

温度 场 的 影响 。 从 图 4 (a) 可 以 看 出 孔径 大 小 对 固 

体温 度 和 流体 温度 都 有 显著 影响 。 孔 径 减 小 时 ， 流 

体温 度 和 固体 都 降低 ， 且 流体 和 固体 间 的 温差 明显 

减 小 。 由 图 4 (b) 可 知 , 当 孔 隙 率 从 0.95 降低 到 0.85 
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图 5 泡沫 参数 对 压 降 的 影响 (a. 孔 径 ，b. 孔 院 率 ) 


Fig.5 Effects of foam structure parameters on 


pressure drop(a. pore size, b. porosity) 

4.4 换 热 特性 

如 图 6 所 示 ， 为 金属 泡沫 结构 参数 对 换 热 系数 
的 影响 。 由 图 6 (a) 可 知 , 管内 充分 发 展 的 Nu 数 随 
着 孔径 的 增 大 逐渐 减 小 。 结 合 图 5 (a) 不 难 发 现在 
在 本 文 的 研究 工 况 下 ， 和 孔径 的 最 佳 取 值 范围 为 
lmm-2mm。 这 是 由 于 当 了 和 孔径 小 于 lmm 时 ， 压 降 急 
剧 增 大 ， 但 是 Nu 数 却 仍 以 相同 的 梯度 增加 。 由 图 6 
(b) 可 知 ，Nu 数 随 着 孔隙 率 的 增 大 而 减 小 。 这 是 由 
于 孔隙 率 增 大 时 ， 有 具有 高 导热 系数 的 固体 相 体积 4 
数 减 少 ， 使 得 有 效 导热 系数 减 小 ， 从 而 导致 换 热 减 
弱 。 
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图 6 泡沫 参数 对 换 热 系 数 的 影响 (a. 孔 径 ，b. 孔 际 率 ) 
Fig.5 Effects of foam structure parameters on heat 
transfer coefficient(a. pore size, b. porosity) 
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5 结论 


本 文 对 完全 填充 金属 泡沫 的 水 平 圆 管 内 的 流 场 、 
温度 场 进行 了 数值 分 析 ， 讨 论 了 金属 泡沫 的 孔径 、 
渗透 率 及 流动 参数 对 流动 和 换 热 的 影响 。 主 要 结论 
如 下 : 

(1) 填充 金属 泡沫 后 管内 充分 发 展 的 速度 分 布 
不 再 符合 抛物 线 型 分 布 ; 管内 流 场 变 得 比较 均匀 ， 
金属 泡沫 圆 管 壁面 处 的 速度 梯度 远 高 于 光 管 ， 壁 面 
边界 层 明 显 变 薄 ; 

(2) 随 着 金属 泡沫 的 孔径 增 大 ， 管 内 轴 向 速度 
增 大 ， 而 壁面 边界 层 厚 度 也 增 大 ; 

(3) 孔径 减 小 时 流体 和 固体 间 的 温差 明显 减 小 
而 孔隙 率 减 小 时 ， 对 流体 和 固体 间 的 温差 影响 不 明 


三 | 
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(4) 孔径 越 大 ， 孔 际 率 越 高 时 ， 压 降 越 小 。 孔 
径 对 压 降 的 影响 大 于 孔隙 率 的 影响 ; 

(5) 孔径 越 小 ， 孔 际 率 越 低 ， 换 热 越 强 ， 但 是 
此 时 的 压 降 也 越 大 ; 

(6) 在 本 文 的 研究 工 况 下 ， 最 佳 的 孔径 范围 是 
1~2 mm。 
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